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The change in spanwiae load distribution at high Mach numbers 
is of considerable interest beoause lateral movement of the center of 
lift affects the trim, stability, and structural factors of safety 
of an airplane. The effects of Mach number upon the span load distri- 
bution of various airplanes have been predicted from the results of 
several series of model tests in high-speed wind tunnels, but all tests 
were not made under comparable conditions.' Different model support 
struts were used and various methods of tunnel calibration were 
employed. ’Fartheimore, the results of the various tests were not 
presented in the same (Banner. This report is an attempt to correlate 
these data by reducing them to similar test conditions and parameters. 

Although the available data are too meager to permit isolation 
of the effects of changes in wing configuration, in a majority of 
cases the tests reveal a tendency for the center of lift to shift out- 
board with increasing Mach number. The exceptions noted are a wing 
with a thiokened and swept-back root section which housed an air-inlet 
duot, and a highly tapered wing of constant thickness ratio with about 
2° of washout. 

V 

The theoretical span load distribution was computed and oampared 
with the wind-tunnel results for the two oases where high-speed air- 
foil section data were available and the agreement was quite, satis— 
"factory# 
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IfflTRODUCTTON 

Most high-speed airplanes attain velocities even in shallow 
dives great enough to cause the critical Mach number of the wing 
to be exceeded. As the critical Mach number usually varies 
across the span of the wing due to changes in wing profile and 
thickness, or to the presence of nacelles, local sweepback, and 
fuselage interference, the load distribution may be strongly 
affected by ccaapressibility. If -the center of lift shifts outboard, 
the bending moment in the wing root may become large enough to 
cause structural failure. Although several series of high-speed 
span— load investigations have been made during recent years, the 
results are not readily comparable because. the models represented 
widely different types of airplanes and were tested on various 
support systems. All basic data in this report are from tests 
in the Ames l6— foot high-speed wind tunnel except those of the 
high-aspect— ratio wing which are from tests in the Langley 8-foot 
high-speed wind tunnel. As improved methods of calibrating wind 
tunnels and - computing cons briction corrections were developed, 
they were incorporated into the test procedures,. Finally, the 
results Of the teats were presented in various ways so direct 
comparisons bf the data were difficult. Since this, report is 
merely a summary of 'the results of tests of the unrelated models 
previously mentioned and includes no further investigations in 
which the effects of the variables were systematically determined, 
most rof the deficiencies of the original results still exist. 
However,- all data have been reduced to, the same parameters and 
the same constriction and calibration factors have been applied 
where possible. 


MODEL AND APPARATUS 

The models tested in the Ames 16-foot wind tunnel were three- 
dimensional models of complete airplanes and were mounted on two 
front struts and a tail strut. A 'sketch showing plan and front 
views of the left half of each model, the location of the front 
support struts, and a table of basic data are included in figuros 1 
through 7. The front support struts for the model of figure 6 were 
so arranged that the leading edges of the struts were at the trailing 
edge of the wing. These struts had a thickness— to-chord ratio of 
0.12 and most of each strut was covered with a fairing. For the 
other Ames tests the front support struts had a thickness— to-chord 
ratio of 0.05 and chords of 8 inches at the model and 6o inches at 
the tunnel wall. The roar strut had a thickness— to-chord ratio of 
0.07 and a 20— inch chord. All three struts were used without fairings. 
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Ble high-aspect— ratio wing modal (fig, l) tested in the Langley 
8-foot high-speed wind tunnel was mounted on a thin vertical 
plate on the model and tunnel center line. 

The locations of rows of wing pressure orifices are 
indicated "by dashed lines on the sketches of the models. 

SEDUCTION AND ANALYSIS GS DATA. 

Symbols 

Bie symbols used, in this report are defined as follows: 

A aspect ratio (b 2 /s) 

a speed of sound in free stream, feet per second 
b Bpan of wing, feet • 

Cp, lift coefficient (L/qS) 

c section chord of model, feet 

7T mean model chord (s/b), feet 
c n section normal— force coefficient (n/qc) 

g acceleration equal to gravity, 32.2 feet per second per second 
L-' lift on model, pounds 

M Mach number (7 /a) 

Mqj. critical Mach number (free— stream Mach number at which the 
local velocity becomes sonic) 

n section normal force, pounds /t=T. 

q dynamic pressure in free stream, pounds per square foot (J-pV 2 ) 
S area of wing, square feet 
V velocity of free stream, .feet per second 

distance along semispan from center line of symmetry, feet 
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a angle of attack of ving chord line or fuselage reference line 
(see figs, 1 to 7), corrected for tunnelrrwall effects and 
upflcw, degress 

p mass density in free stream, slugs per cubic foot 


Methods of Reduction 

In every case the pressure-distribution data were reduced 
to section normal— force coefficients by integration of ourves 
giving the upper and lower surface pressure coefficients along 
the chord of each section. No corrections for wing dihedral were 
made to the section loading coefficients c'nc/c’ when they were 
plotted along the model span. The centroids of the loading were 
determined by integration of curves of section loading coeffi- 
cients plotted along the span. 

The lateral oenters of lift were determined by integrating the 
loading curves shown in figures 1 to 7 without extrapolating the 
curves inboard of the points at which the data terminate. Therefore, 
the lateral centers of lift shown are actually the centroids of the 
loading outboard of the pressure stations at which the loading 
curves terminate and show the effects of Mach number upon the bending 
moment at these stations rather than at the actual wing root. 


Methods of Analysis 

In this report the theoretical loading was computed by a 
generalized method of applying lifting— line theory (reference l) 
which utilizes actual high-speed wind-tunnel data. Most previous 
reports on span-load tests presented a theoretical loading computed 
by the method of ANC— 1(1) (reference 2) using experimental or 
estimated low-speed airfoil characteristics. 


RESULTS AND DISCUSSION 

The Ames l6-foot wind— tunnel teste summarized in this report 
cover a Mach number range of 0,40 to 0.80 or 0.825 and the Langley 
8— foot wind-tunnel tests extend from a Mach number of 0. 40 to 0.925* 
The Reynolds number for the Ames tests varies from 4, 000, 000 for 
the smallest model at 0.40 Mach number to 9,000,000 for the largest 
model at 0.80 Mach number. Reynolds numbers of 900,000 at a Mach 
number of 0.4 and 1,400,000 at 0.90 Mach number were attained during 
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the Langley tests. 

The variation of critical Mach number along the model span is 
shewn in figure 8 for several of the models tested. Models A and C 
shew the normal variation, the thicker root sections having the 
lower critical Mach numbers. Quite the reverse is true in the case 
of model B where the inboard section has the highest critical Mach 
number except at negative or large positive -lift coefficients when 
the pressures reach critical values over the duct lip, This unique 
effect appears to be the result of several. factors. When air is 
flowing through the duct a higher critical Mach number is maintained 
over a wider range of lift coefficients than when' the ducts are 
closed. Even with the ducts closed, however, the root section has 
a higher critical Mach number than an HACA 66-series airfoil of 
comparable thickness— to-chord ratio; so it is obvious that the 
excellent high-speed characteristics are due in pert to three- 
dimensional effects produced by the marked sweepback and taper 
existing over this portion of the wing. Model D has a relatively 
small variation in critical Mach number across the wing span, the 
root and tip sections having somewhat higher values for all lift 
coefficients up to about 0.70. A 3 csnewhat similar variation of . 
critical Mich number exists across the span of model E. Hie tip 
sections always have the highest critical Mach number and the 
middle portions of the wing have the lowest. The slight Increase 
in critical Mach' number at the point of dihedral reversal 
(approximately Q.k semispan) for lift coefficients . of 0.60 and 
0.80 is probably due to local separation and the reduced lift 
coefficient occurring at this point. 


Lift Characteristics 

The variation of the lift— curve slopes and angles of zero 
lift with Mach number is shewn in figure 9# In every case the 
lift— curve slope increases smoothly up to some Mach number between 
0,675 and 0.775 and then breaks sharply downward. This Mach 
number of divergence varies with the wing thickness, the break in 
the curve occurring at the lower Mach numbers for Hie thicker wings. 
The angle of attack for zero lift remained constant or decreased 
slightly up to about 0.70 Mach number for all models. With further 
increases in Mach number the zero— lift angle -of models B, D, and E 
shifted abruptly in a positive direction, while the increases for 
models C and F were more gradual. Up to 0.80 Mach number the 
zero— lift angles of Model A and the high-aspect— ratio wing had 
changed only about 0.5° although the high-espect— ratio wing 
experienced a rather large increase above 0.825 Mach number and 
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it appeared that model A might do likewise at the higher Mach ' 
numbers. Model B experienced . the largest increase in zero— lift 
angle, the shift being nearly 4° between Mach numbers of 0,40 
and 0 . 80 . 


Span Load Distribution 

As the models for which high-speed span— load data are 
available vary widely in general configuration, airfoil section, 
thickness and twist distribution, plan form, and aspect ratio, 
it is not possible to isolate the effects of these variables. 
Theoretical analysis of the loading is further complicated by 
the fact that most models employed modified airfoil sections for 
which no high-speed section data were available, Consequently, 
all that could be done was to reduce the original data to similar 
test conditions and to the same loading coefficients and to 
present it at the same values of angle of attack and Mach number. 

Although the low— speed span loading can usually be predicted 
with reasonable accuracy using the method of ANC— l(l) (reference -2), 
this procedure is of little value for high-speed conditions 
because the lift characteristics often change considerably and 
even become nonlinear. Therefore, the loading for the high-aspect— 
ratio wing was computed hy the method outlined in reference 1 , 
using actual high-speed section data, and the results are 
presented for 0.40 and 0.80 Mach number in figure l(a). This wing 
was particularly amonable to theoretical analysis since it had a 
simple plan form and an undistorted airfoil section for which 
high-speed section data were available. It is apparent that the 
theoretical and measured load distributions are in excellent 
agreement when compared at equal lift coefficients. The loading 
changes quite radically at Mach numbers of O .85 and above. From 
figure 10 it may be seen that for all positive lift coefficients, 
the center of load of this model moves outboard above 0.77 Mach 
number. The shift is greatest at low angles of attack and would 
result in variation in trim and stability due to changes in the 
downwash at the tail,- A further effect of this shift in the 
center of load is to increase the root bending moment . For a 3s 
pull-out the increase is 5 percent as the Mach number changes 
from 6,77 "to 0,90. This banding moment is about 2.5 percent 
greater than that predicted by ANC— 1 ( 1 ). 

Hie loading on the wing of model A, shown on figure 2, is 
quite different from that on the high— aspect-ratio wing. Part of 
this difference is apparently due to a substantial amount of lift 
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contributed by the fuselage because the loading at the wing root 
is too great to be accounted for by the additional area or the 
change in effective, angle of attack caused by the small leading- 
edge extension at the wing root. Although there is a srmal 1 
increase in root bending moment due to the outward shift of the 
center of lift with increasing Mach number, the root bending 
moment never quite attains the value predicted by ANC— 1(1). Since 
no high-speed section data were available for the airfoil used on 
this wing, no other theoretical calculations were attempted. 

Model B (fig. 3) is another example of an airplane with a 
wing which departs radically from a straightforward design. In 
this case however, it is the unusual wing root which carries more 
than the expected load at high Mach numbers. The experimental 
results indicate tip loads lower than. those computed by ANC— 1(1) 
and for this reason the root bending moments - are less than the 
theoretically predicted oneB at nearly all lift coefficients and 
Mach numbers. The somewhat . irregular distributions measured at 
the higher Mach numbers may be due go support strut interference. 

Model C (fig. 4) has a wing of conventional plan farm, 
thickness variation, and twist distribution. With increasing 
Mach number the lift for a given angle of attack first increases 
as a result of the increase in lift-curve slope, Due to its 
greater fchickneas and higher angle of attack the root section 
exceeds the critical Mach number first, resulting in a decrease 
in load carried- by this section. At the higher Mach numbers sill 
sections of the wing are operating above their critical conditions 
resulting in a general decrease in loading across the span. The 
loss of lift at the. root section causes an outward shift in the 
center of loading with a consequent increase in wing root bending 
moment. At a Mach number of 0.80 and a lift coefficient corre- 
sponding to a 3g pull-out at 5000 feet the bending-moment at the 
wing root is about 17 percent greater than for a similar pull-out 
at 0.6 Mach number. 

Although the wing of model D (fig. 5) and the high-aspect— 
ratio wing (fig. 1) are quite similar in plan form they differ 
considerably in sweepback, dihedral, . thickness, and twist. Also, 
the wing of model D was tested in combination with a fuselage. 

The thickness and twist distribution of the model D wing appear 
quite satisfactory, since the wing does not experience any large 
loss in lift over the root section at the higher Mach numbers. 
Figure 10 shows' that the lateral center of lift actually shifts 
inboard with increasing speod. ' 
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Model E (fig. 6) is a special case, since it is the only model 
of the group with a wing of the inverted-gull type. The large loss 
of lift over the inboard section of the wing is due to the thicker 
inboard wing panel having tfce lowest critical speed and to separa- 
tion of flow at the point of dihedral reversal,.. Figure 10 shows a 
continuous- outward shift of the lateral center of lift as the Mach 
number is increased. This shift of loading results in an increase 
in bending' moment at the wing root of approximately 15 percent 
above the theoretical value for a 3s pull— out at a Mach number 
of 0.8o and an altitude of 5000 feet, 

Although' rather incomplete span-load data were obtained 
during the model F tests, because the model had but three wing 
pressure stations, the data have been included because the simple 
and straightforward design of the wing make it- amenable to 
theoretical treatment. Whiie the airfoil employed is not a basic 
section, it is closely related to one for which high-speed section 
data are available. The upper sets of curves in figure 7 show the 
span loading at 0.40 and 0,80 . Mach number as . determined from wind- 
tunnel tests. The lower sets of curves are for the loading chaaputed 
from section data by the method of reference 1, and the test points 
shown are from measured loadings at the same lift coefficient as 
the computed loadings. As no wind-tunnel data were available for 
the 65 (112) —21 3 airfoil the section characteristics used in the 
computations were obtained by plotting available wind-tunnel data 
for the 65]_— series airfoils and then extrapolating to got 
characteristics for a 651— 213 airfoil. The agreement between the 
theoretical and measured loadings is excellent at 0.40 Mach number 
but the loss of lift measured at the wing root is not indicated by 
the theoretical analysis. It is believed that this loss of lift 
at the root may be the result of interference between the wing and 
fuselage because the fuselage air— intake- ducts undoubtedly have an 
appreciable effect upon the wing— root pressures. 


CONCLUDING REMAKES 

Correlation of available span— load data from high-speed wind- 
tunnel tests indicated that in the majority of cases the center 
of lift shifted outboard at the higher Mach numbers. Although this 
shift resulted in wing— root bending moments greater than those 
indicated by low-speed wind-tunnel tests, in some cases the center 
of lift did not shift outboard of the point computed by low— speed 
theory. For two of the models tested the center of lift actually 
shifted inboard with increasing Mach number. The excellent high- 
speed characteristics of the one model, which had a thickened wing 
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root housing an air-inlet duct, were apparently due to three— 
dimensional effects produced by the sweephack and taper of the 
root section. The other model for which the center of lift shifted 
inboard had a highly tapered wing of constant thickness ratio with 
about 2° of washout. The outboard shift of loading for a model 
with an inverted-gull wing was apparently due to prematura separa- 
tion and loss of lift at the point of dihedral reversal. 

When high-speed airfoil data were available and the model wing 
was of conventional plan form and was not subjected to large ■ 
interference effects from other parte of the model it was possible 
to compute the span load distribution at high Mach numbers with 
considerable accuracy. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif. 



MCA EM No. A7C28 


EEFISENCES 

Boahar, John: The Betermination of Span Load Biatrihution 

at High SpeedA-hy IJae of High-Speed Wind-Tunnel Section 
Bata. NACA ACR 4B22, 1944. 

Anon.:' Spanwiee Air-Load Biatrihution, ANC— l(l). Army— 
Navy-Commerce Committee on Aircraft Requirements. 

April 1938. 



Fig. la 




f<r) M • O.4oo 7b M r o.&oc 


F/avRG A - Spam ms e vak/a-tto*/ or sscr/o*/ loao/ms for a 
M oose, or a Hrcsrt • A Spec t~- ra no ma/«. 


NATIONAL ADVISORY 
COMMITTEE POt AERONAUTICS 























Fig, lb 


NACA RM No. A7C28 




o £ * Ji A AO 

£y/b 


o 2 * .6 .a ao 

2y/h 



JvGu/ZE A- Couccu£>eG . 


NATIONAL ADVISORY 
COMMITTEE FOR AERONAUTICS 


NACA RM No. A7C28 


Fig. 2 



NATIONAL ADVISORY 
COMMITTEE FOE AEBOWUTICS 
































Fig. 3 


NACA EM No. A7C28 



NATIONAL ADVISORY 
COMMITTEE FOR AERONAUTICS 













NACA RM No. A7C28 
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